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I. Introduction 


The advent of high speed digital computers has made a striking im- 
pact in aircraft flight control systems. In the flight control actu- 
ation systems alone, digital technology has brought about the possi- 
bility of quite a number of new design approaches [1]. For instance, 
the digital actuator is one of these new possibilities. Basically, a 
digital actuator is a device that can produce a desired quantized rate 
or a quantized position. The question that immediately arises is 
how often the on-board control computer must output digitized rate or 
position commands to the control actuators. If a fixed high output 
rate is selected, the control system analysis is greatly simplified 
enabling the designer to use the standard discrete control theory. On 
the other hand, a fixed high output rate increases the burden on the 
flight computer and the wear and tear of the associated mechanical 
equipment. Intuitively, a multirate output sampling is needed to im- 
prove the sampling efficiency. In a multirate sampled digital environ- 
ment, the control computer outputs per unit time can be optimized to 
achieve a desired system performance with the least number of output 
commands . 

The desire to achieve a high level of sampling efficiency in sam- 
pled-data control systems has been a challenging goal in control engi- 
neering. Many different methods have been used to obtain efficient 
adaptive sampling laws in the last decade [2] - [10]. In most of these 
methods, an integral difference criterion is used to make the difference 
of any two consecutive sampled errors equal. In general, these methods 
have been relatively simple formulations and yield implementable laws 


with low computational requirements. The major drawback in all of these 
approaches has been the lack of guarantee of stability in the resulting 
closed-loop system. One reason for this stability problem is that, in 
all these systems, the sampling interval is selected from a continuum 
of sampling intervals which makes the problem nonlinear. In this study, 
the sampling interval is selected from a finite number of sampling in- 
tervals . The use of a finite number of sampling Intervals results in 
an analytically tractable formulation and also fits better to physical 
reality since the discretization of the continuum of sampling intervals 
is necessary in order to be implemented in a digital environment. 

The outline of the report is as follows: In the next section, a 
mathematical statement of the problem is given. Basic notions of sta- 
bility for time-varying systems is reviewed in Section III . In Section 
IV, it is shown that the stability problem as modelled in this study es- 
sentially reduces to one of finding a convergent sequence of matrices 
from a finite collection of matrices. The notions of convergent, con- 
tractive, and pre-contractive set of a finite number of matrices is intro- 
duced. It is shown that pre-contractiveness is a necessary and sufficient 
condition for finding a convergent sequence from a finite set of matrices. 
In Section V, computable criteria are found to test the contractiveness 
of a set of matrices. Section VI contains the conclusions and directions 
for further research. 
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II . S tatement of the Problem 


The objective of this investigation is to find mathematical con- 
ditions which will ensure stability in linear multirate-sampled-data 
control systems. Specifically, we take as given an m x m constant 
matrix A, and an m x n constant matrix B, and consider the system 


Cl) 


X = Ax + Bu 


where x is an m-component "state vector" depending on time, and u is 
an n-component "control vector" depending on time. Given an initial 
state Xq = x(0) , we wish to construct a control u which will produce- 
stability of the system at zero. (Various types of stability will be 
discussed in the next section.) 

All numbers, vectors, and matrices considered in this report are 
assumed real. Complex numbers enter only in the brief allusions to 
the spectral radius of a matrix in section 4. 

We must first discuss a solution of (1) satisfying an arbitrary 
condition of the form x(t]^) = z, where tj^ > 0 and z e r”'. The fol- 
lowing solution to this problem is standard. 

Consider the equation x = Ax. Any m x m matrix $ whose columns 

are linearly independent solutions of x = Ax is called a fundamental 

t A. 

matrix for the equation. One such matrix is $(t) = e . Furthermore, 
if t^ > 0 and z e R™, a solution <j) of x = Ax satisfying if>(t^) = z is 


given by 


<),(t) = e^"-^l^^ 
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Also, a solution 9 of the original equation (1) satisfying 0(t^) = 0 

is given by ft . 

0(t) = $(ty $”^(s)Bu(s)ds = e 1 e ®%u(s)ds . 

Jti 

Now (ij) + 0) ' Ct) = A(j)Ct) + A9(t) + Bu(t) = A((j) + 0)(t) + Bu(t) and 
(<)> + 9)Ct^) = <|>Ct]^) + 0Ctj^) = z. Thus X = (j) + 6 is a solution of (1) 


satisfying "xCt^5 = z. 

Now letting T = t-t^^, we may write x(t) = e^^ ^ 
e^^z + e“^®'''*^l^^uCs+tj^)ds = e^^z +J^ e^'^"®^^u(s+t^) ds . 


^)A tA/ ( 
1 z + e J 


: ®^u(s)ds 


Substituting v = T-s, we obtain 

(2) x(T) = e"^^z + y* e'^^u(t^+T-v)dv . 


We now think of constructing a "stable" solution of (1) by the 
method of sampling; that is, at particular instants in time we examine 
the value of x and, on the basis of that value, change the value of u 
so as to make x stable. The value of u is to be held constant until 
the next sampling instant. Thus if Increasing sequence 

of sampling instants, Tj^ = - tj^ is the k*^’^ sampling interval, 


Xj^ = xCtj^I is the state at time tj^. 


and 




is the constant value of 


u(t) for t, < t < t, , - , formula (2) can be rewritten as 
' ' k k+1 


^+1 


TkA , 

= e k X, + 




fTK . 

- 


\ 


dv 


or 


(3) 


X, = e^^^ X, + 

^+1 * \ 


i 


vA 

e dvB u^ 


A variable rate sampling problem is a problem of the above type 
in which a collection S = -js^, S^, • * *, of positive numbers 
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(sampling interval lengths) is given. The sequence sampling 

instants is generated in time along with the sequence of controls j 

in that each time a sample Xj^ is observed, we are allowed to use its 
value in determining the next sampling interval, Tj^, selected from the 


set S . 


The next control, Uj^, is to be determined in the following manner. 
Corresponding to each in S we determine an n x m matrix called 
the i^^ feedback matrix . Whenever the sampling interval T^^ is selected 
to be the member of S, then u^^ is determined by u^^ = Further, 

we wish to allow the choice of T^^ to depend not only on the observed 
state But also on the previous sampling interval 

Thus, given A, B, and S, we wish to define the feedback matrices 


^2’" ' '» ^ "decision function" ^k-1^ whose value is 

a member of S which we will choose as T, . In this way, starting from 

tC 

any initial state x, and an arbitrary member of S, a sequence of 
state values is defined, using (3) by: 



We want to choose the F^'s and the d in such a way that the resulting 
system is stable at 0. 
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Ill . Stability 


We restrict our attention to notions of stability for sequences 
of vectors in R™ generated by some discrete system [11]; that is, 

^+1 = » ^1 

In general, P may involve any number of parameters such as, in our 
problem, The notions of stability require a norm on R^. It 

is the case that all norms on R™ are equivalent in the sense that , 
if II • 11^ and | | • 11^ are any two norms, there are numbers a > 0 
and S > 0 such that 

M 5K X Il 2 and || x ||2<3|| x ||^ 

for all X e k™. As a result of this, it is immaterial which norm is 
used in the following definitions, and we denote a general norm by 

II • II. 

Several definitions for various types of stability have been 
defined. The ones given here essentially follow Ortega [12], 

The discrete system P is convergent of x* if there is a 6 >0 
such that 

lim I I X, - X* I I =0 whenever | | x^ - x | | < 5 . 

■k 

P is globally convergent at x if 

II ^11 tn 

lim I I X, - X I I = 0 for all- x e R . 

•k 

P is Lyapunov stable at x if, given e > 0, there is a 6' > 0 
such that I I Xj^ - x | | < e fo£ all k > 1 whenever | |x^ - x | | <6. 
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Global convergence does not imply Lyapunov stability, for 
1 * 

consider In R , with x *• 0, 

\+l “ “ k ’ ^ ^ ^1 ^^^itrary. 

Also, Lyapunov stability does not Imply convergence, for consider In 

1 * 

R with X •• 0, 

! 2 x ^, k=*l 
x^, k>l • 

* 

P is asymptotically stable at x if P is both Lyapunov stable and 

A 

convergent at x • 

P Is expontlally convergent at x If there is a T >0, a positive 
g < 1, and a 6 > 0, such that 

I |xj^ - X II < rg^l |x^ - X II for all k > 1, whenever 

I I ’'l - X* I I 5 . 

■k 

Global asymptotic stability at x is Lyapunov stability together 

* 

with global convergence at x , and global exponential convergence is 
defined by requiring 

I |xj^ - X II < rg^l |x^ - X I I , all k s 1 and all x. 

Clearly (global) exponential converge implies (global) asymptotic 
stability . 
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IV. Converaent, Contractive, and Pre-contraetive Sets . 


In order to attack the problem as stated at the end of section 2, 
we proceed as follows. Assuming the m x m matrix A, the m x n matrix B, 
and the set S = S 2 , * • *, of positive sampling Interval lengths 

are given, set 


r* ^ i^ j 

C = e , and 


"i-f 


vA 

e dvB . 


If Efc is to be the feedback matrix corresponding to S . , equations (4) 
^ 1 

at the end of section 2 can be written: 




- (S ^ • 

If suitable feedback matrices and a suitable decision function are to 
be found, then surely the collection 


"1 


H. = C. + D.F , 1 < i < 

X 1 IS. 

X 


j 


has the property that, given any vector x e R , there is an order in 
which we can apply the to x to produce a sequence convergent to, 

or at least remaining near, zero. We will find conditions which assure 
that a set of matrices has this property, in the hope that then the F 

°i 

can be found to produce these conditions. 

In our setting it will be clear that convergence, exponential con- 
verge, and asymptotic convergence are equivalent to the corresponding 
global types of stability. In the remainder of this report, the adjec- 
tive "global” will be dropped. Also, all types of stability will be 
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if if . 

considered at x = 0 , so the phrase "at x " will be dropped. Through- 
out the report, theorems are numbered consecutively , regardless of the 
section in which they appear, as are definitions, lemmas, corollaries, 
and numbered formulas. 


Definition 1 : The set K = ’ • ’ , 


of m X m matrices 


is convergent provided that, if x e R , there is a sequence 
with p(x)^ e |l, 2, • • •, n} for all 1, such that 


{p(x)i}“^l, 


lim "tTH z' ^ .= 0 . 


(The product is taken in reverse order so that , is applied first, 

H , ^ second, and so on.) 

P Cx) 2 

In the familiar case where N = 1, the only possible sequence is 


|«l^-|k=l 


. It is well-known that this sequence converges to zero for 


all X if and only if the spectral redius of is less than 1, where 
the spectral radius is defined as the largest number which is the ab- 
solute value of an eigenvalue (real or complex) . That this does not 
generalize in any foreseeable way to N > 1 is seen by the following 
example in which N = 4: 


Example I: 


1/2 0 


1.25 -.75 

-.75 1.25 


0 1/2 


1.25 
.75 1.: 


In this example, each has spectral radius 2, and yet the set 
K = is convergent. This can be seen by observing that: 

If I 1 • II denotes the Euclidean norm, | |h^x| | ^ 3| |x| | for all x 0 
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To prove this, let a = — — and consider the function 
f(t) = (l+t2) - (1/4 + 4t2) = 3/4 - 3t2. f(0) = 3/4, f(l/2) = 0, 

f decreases on ^0, 1/2], and f is even. Since 0 < tt/ 8 < 1/2, 
f(Tr/8) > 0. fCir/S) = 3/4 - 3ir^/64 = a. Assume |t| i tt/8. 

Then f(t) S: fCTr/8) = a. Thus (l+t^) - (1/4 + 4t^) > a 


, 1/4 + 4t^ 

^ i ' + 


1+t^ 


1/4 + 4t^ 

—TTV~ 


< 1 - 


i+t 


2 - < 1 - 


l+ir^/64 


1 _ 

64+ir^ 


(5) 


• • • < e for Itl < I 




+ t^ 


Now suppose Cx,y) e and ^ ^ • 

iX O 
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■ 1 * 1 ^ + 

llgjil - u\f^ 

Thus by- QS) , 


|Hl 


II 

II 


II 


Ji+ 4 (i)^ 

V /fr ^ ®- 

^ (f 


and 


Since C- 


l»i^ll * e||[^l| . 

{(,.,)} I li^l . I . Hj, 


the statement Is justified. 


Furthermore H^, and contract the lengths of vectors in C^, C^, 

and in the same way, where is the counterclockwise rotation of 
through 45°, is the counterclockwise rotation of through 90°, 

and is tEe counterclockwise rotation of through 135°. 

AO 2 

Now = R . Thus, given x e R , we select p(x)j^ = j if 

have been determined, we select 3 

(If at some stage there are two possible choices 


X e C. and, if 
. 1 ' 

if. 




for j, that is two cones overlap, we may choose either.) 


Then 



< e^'x 


0 as k — ^ ' 


This example suggests the following definition. 

Definition 2; A set K = • • •, of m x m matrices is contractive 

relative tc tne norm [ | • | | , ^ provided that, if x e R™ and x ^ 0, there 
is an i e ^1, 2, •••, N j such that | 1 h^x| < ||xl|^. 

It is a theorem that if K is contractive relative to some norm, 
then K is convergent. This will be contained in the more general 
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Theorem 1, to Be proved presently. 

Clearly, contractiveness relative to a given norm is not a 
necessary condition for convergence of a set K of matrices. For 
example. 


K: 






0 ol . 

jo 3J ^ 


xs not contractive 


relative to the Euclidean norm, for neither matrix contracts 
X = Cl,l). Yet the matrix is zero, so that, if 





i=l 

i>l 


for all X, 


then all x go to zero in a hurry. In fact, it is clear that, if 

and some finite product of the H^'s has spectral 
radius less than 1, then the sequence p(x) can be chosen to produce 
that product repeatedly, and thus K will be convergent. 

The authors have not proved that no finite product of the four 
matrices in Example 1 has spectral radius less than 1, but extensive 
experimentation suggests that this is the case. Yet the set K of 
Example 1 is convergent, and the algorithm for selecting p(x) does 
not deal with the spectral radii of products, but rather with the 
more near-sighted contraction of the length of the vector. 

Consider another example with N = 2. 

Examp le 2 : 


K: 


H 


1 


1/2 0 
0 2 



/3?2 1/2 

- 1/2 


is identical with the of Example 1, and effects rotation 
o 2 

clockwise through 30"'. Clearly K is convergent, for, given x e R , 
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we apply H 2 often enough to bring the vector into the cone 
described in Example 1, and then apply Since | lH 2 y| | = | |y| | 

for all y, we have now contracted the vector. If the result of 
puts the vector outside of C^, we again apply H 2 often enough to 
bring it back to C^, then apply and so on. This example, to- 

gether with, the foregoing discussion, suggests the following defini- 
tion. 

Definition 3 : A set K = ^ of m x m matrices is pre-con- 

tractive relative to the norm [1-1^ provided that, if x e R™, x 0 , 
there is a finite sequence s |l, 2 , • • ■ ,n|- , such that 

l(fr "qCrij- 

i=n(x) 



We now show that pre-contractiveness relative to some norm is equiva- 
lent to convergence. In the following discussion, K = Jh , H X 


is a set of m x m matrices and 


■ V 


is a norm on R . 


Lemma 1: If K is pre-eontractive relative to | | • | i , there is a posi- 

V 

tive integer M such that, if x e R*^, x 5 ^ 0, then there is a finite 
sequence e |i,2,*'*,n|, such that 


K’-fr 

i=n(x) 


X < 

V 


' V 


and n(x) < M. 


|x e r” 1 |xl 1 ^ = l|. 


Proof: Let B = £ R 

B Is compact in the norm topology. For each positive integer n, let 




S = /X e B 3 q(x) 


®q(x)^V 

f=k ' 


< 1 , and k 
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S is open relative to B and u,S = B, since K is pre-contractive 
n n=l n ^ 

relative to |1*||^. Thus there are positive integers n^^, n^, •••, n^ 


such that 


US = B. 
i=l “i 


Let M = max D 2 ’ **’’ 


Given 0 * x e R™, there isanj,l^j^)l such that T^-~r~ ^ S 

I • 

Then 


'k'im) 




V 


|x|l^ • 1 = ||x| 1^ and n < M . 


QED 

I I , there is a 3 , 


i e -|l,2, • • • ,N 


Lemma 2 : If K is contractive relative to |{. 

0 < 6 < 1, such that, ifxeR”*, xa^O, then there is an i 
such that I |h^x| |^ < 3||x| |^ . 

Proof: Let B be defined as in Lemma 1. For each positive integer n, 


let 




X e B 


3i e 




oo 

Each S is open relative to B, and u-S = B. Thus there are positive 
n n=l n 

integers n^, n 2 , n^ such that 


,u^S = B. 
i=l ni 


Let n« »= max /n, ,n„ , • • • ,n. I, and let 3 = — . Then 0 < 3 < 1. and. 

0 \ 1 2 Zf nQ » 


HI X 

given X e R , x ^ 0, choose j so that 


, . I I e S , i.e. 


H, 
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Then 1 |Hjx| 




|x|l^'’ - ^11^1 


'V 

QED. 


Theorem 1 : K is convergent If and only if K is pre-contractive rela- 
tive tq 1 I . 1 ly. 

Proof ; Assume K is convergent. Choose x e r’”, x 0, and let p(x) 
be a sequence such that 




lim , 
k->«\i 

There is an n(x) such that 






I X I I ^ , and we may 


choose q(x)^ = p(x)^ for 1 < i < n(x) . Thus K is pre-contractive 
relative to I I • I I . 


Now assume K pre-contractive relative to | | • | |^. 
positive integer M as given in Leirana 1. Then the set 


Choose a 


/V j 1 

K = »TTH 


) 1 1 “ 

(i4 Pi 

is contractive relative to 


k < M 


, p^ e |i,2,"‘,n||' 


. Choose B, 0 < 3 < 1, for K 


as 


given in Lemma 2 . Thus we have : 


If X E R™, X 0, 3|q(x) ^ 

I 1=1 

.11 (it ^%Cx) .)^1 Iv “ ®1 l^'l U and n(x) < M. 


i=nCx) 


Let 


be the operator norm on m x m matrices generated by | | ' | |^> 
^1 Ill’ll 1^ Iv " l}- I |Hy| 1^ ^ I |h| IqI |y| 1^ for 


i.e. H L = sup< Hx 


all y E R®, and 1 Ih.H^] |q < j | h^| |q1 \n.\\^. 

{iKIlo, l|H,llo. IIVIo} 


Let B > max 
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Fix an 


X e R™. We must produce a /p(x)l so that lim(7TH . v lx = 0. 

\ k->coW4 


If X = 0, any p(x) will do, so assume x ^ 0. 
Define = x and, for k > 0, 


k 

y = 




k-1 


k-1 k 

unless y = 0, In which case let y =0. Then 

For each positive Integer k, let iDj^ be the unique non— negative Integer 
satisfying : 

i\-l . 

s(k) = IZ n(y-') < k and t(k) = Hn(y^) > k, 

j=0 


j=0 


(where we may take n(y^) = 1 If y'' = 0) . Then > nij^ for k > 0, 

and lim 




k -^0 


™k 


pCx)j^ Ck-1 ' Q, cKoose a positive Integer Si 

so that 3™^ < — 


|x||^ • B 




Suppose k > Then k-s(k) s t(k) - s(k) = n(y™k) < M. Thus 

."kl 


(ii 

(a") 






m ,. 


,.“k| 




TTB B = 6 '^MxM B“ < e 

i=M ' 


Therefore lim 
k' 


~ ^ ^ convergent. 


QED, 
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We list three corollaries to this theorem. The first involves 


the concept of exponential convergence, which, in this setting, is 
defined as follows: 

Definition 4 : A set K=|h^,H 2 , ' * • ,H^ | of m x m matrices is exponen- 
tially convergent provided that there is a norm | | • | |^ on R™, a 
r > 0 and a 6, 0 < g-< 1, such that, if x e R™, there is a sequence 
p(x)i|i_i, p(x)^ 0 |i,2,‘**,n| such that 






V 


< FB 1 lx| I for each positive integer k. 


Corollary 1 : If K is contractive relative to | | • | |^, then K is expo- 
nentially convergent . 

Proof : Since K is contractive relative to | | • ] | , the M in the proof 

of Theorem 1 can be taken as 1. Thus each n(x) is 1 and it follows 
that nij^ = k-1 for k > 0. Therefore the last inequality in the proof 
becomes : 



< B 


k-1. 


= FB^I |xl 1^, where T 


B 


QED. 


Corollary 2 : If K is pre-contractive relative to | | ' | |^» then K is 
pre-contractive relative to any other norm | ] • | |^ on R™. 

Proof: K pre-contractive relative to |1*||^ implies K is convergent, 
which Implies that K is pre-contractive relative to | ] ’ | |^. 

QED. 


Thus pre-contractlveness is norm-independent and we have: 

Definition 3'^ ; K is pre-contractive if there is a norm | ] * | |^ on R™ 
such that K is pre-contractive relative to | 1 • | |^- 
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Theorem 1 can now be restated: 

Theorem 1’ : K is convergent if and only if K is pre-contractive . 

In the present setting, the notion of asymptotic stability is 
given by: 

Definition 5 : K is asymptotically stable provided that there is a 
norm | | • | |^ on R™ such that, if e >0, there is a 6 > 0 such that 
if X e R™ and ( |x| < 6 then there is a sequence |pCx)|% such 

that 

lim / . . \x = 0 and I 1/ TTH / s \x < e 

k->«U.kP«J ||\i4 '' 

for each positive integer k. 

Corollary 3 : K is convergent if and only if K is asymptotically 
stable. 

Proof : Clearly, asymptotic stability implies convergence. Assume, 
then, that K is convergent. Let | |* | |^ be a norm. K is pre-con- 
tractive relative to | 1 • | |^. Given f > 0, choose M, 3, and B as 
in the proof of Theorem 1. Let 5 = ^ , and suppose x e R™ with 
I I x| 1^ < 6 . Choose ^ as in the proof and we 

have, for each positive Integer k, 

II "fTH x||^ < 6"^l|x||^B^ < = E. 

i=k i 

Therefore K is asymptotically stable. QED. 
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V. Criteria for Contractiveness — Positive Definite Seta of 
Matrices . 

We now seek criteria for contractiveness and pre-contractlveness 
of a finite set of m x m matrices . We have some information con- 
cerning contractiveness for which we must define the notion of a posi- 
tive definite set of matrices . Recall that a symmetric m x m matrix 

T 

A is said to be positive definite provided x Ax > 0 for all non-zero 
X in R™. Equivalently, A is positive definite provided that all the 
eigenvalues of A, which must be real because of the symmetry, are 
positive . 

Definition 6 ; A set W = |Aj^,A 2 , ' • • .A^^^ of symmetric m x m matrices is 
a positive definite set provided that, if x e R*°, x ^ 0, then there is 
an i £ |i,2,***,n| such that x^A^x > 0. 

We will show that a finite set of m x m matrices is contractive 
relative to a certain norm if and only if a certain related set of 
matrices is a positive definite set. 

Definition 7 : If B is a symmetric m x m positive definite matrix and 
X e R™, the B-norm of x is defined by 



It is well-known that | | • | |g is a norm on R™. Notice that B = I gives 
the Euclidean norm. 

Lemma 3 : If H is an m x m matrix, B an m x m positive definite matrix, 

rp 0 T* m 

and A = B - H BH, then 1 |x| [g - | |Hx| jg = x Ax for all x e R , and 
so |x E R™ I I Hx I I g < I I X 1 1 g| = |x e R*” x"^Ax > 0 1 . 

Proof ; Let x e R™. Then | jx| ]g^ - [ jHx| |g^ = 

x^Bx - (Hx)'^BHx = x^(B - H^BH)x = x^Ax. 

QED. 
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Definition 8 : If K = | H^,H 2 , * • * is a set of m x m matrices and 

B is a positive definite matrix, the B-symmetric set of K is the set 


B - H. BH. 

I 1 1 


1 < i 


< n| 


of symmetric m x m matrices . 

In [13], Stein proves: 

Theorem : If H is a real or complex square matrix, a necessary and suf- 
ficient condition that lim = 0 is that there exist a positive 

n-J“ 

definite Hermltian matrix B for which B - H BH is positive definite. 

If H is real, B may be taken real and symmetric. 

Since contractiveness relative to some norm implies convergence, 
the following theorem contains a generalization of one direction of 
Stein's theorem. 

Theorem 2 : If K = |h^ ,H 2 , • • * and B is a positive definite matrix, 
then K is contractive relative to [I'ljg if and only if the B-symme- 
trlc set of K is a positive definite set. 

Proof : Let = B - so that W = ,A^ , • • • , is the B-sym- 

metric set of K. The following statements are equivalent: 

K is contractive relative to B. 

X e R™, X ^ 0, implies | ]h^x| | g < | |x| |g for some i. 

X e R™, X 0, implies x'^A^x > 0 for some 1. 

W is a positive definite set. 

QED. 


Theorem 2 together with Corollary 1 yields : 

Corollary 4 : If K = |h^,H 2 , ’ * * ,H^|, B is a positive definite matrix, 
and the B-symmetric set of K is a positive definite set, then K is 
exponentially convergent. 

Thus we would like to know conditions related to positive definiteness 
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interior 

cone 


of a set of symmetric matrices . 

Recall the following theorem from matrix theory: 

Theorem: If A Is a symmetric m x m matrix, there Is a unique dla- 


|l,l,--,l,-l,-l,..-,-l,0,0,’--,o} 


gonal matrix D of the form D = dlag 

such that there is a non-singular matrix P with D = P'^AP. The matrix 

D can be computed from by several standard techniques, one at which 

is to repeatedly apply similar row and column operations to A. The 

product of the corresponding elementary column matrices gives a 

T 

matrix ,P such that D = P AP . 

Definition 9 : In the preceding theorem, D will be called the canon- 

T 

ical form of A , any non-singular matrix P satisfying D = P AP will 
be called a canonizing matrix for A, the number of I's and -I's on 
the diagonal of D will be called the index and co-index of A, res- 
pectively. 

The positive definiteness of a set of matrices will be studied 
by considering Images under linear transformations of certain fixed 


define 

these sets. 


< p, 0 

^ n, and p + 

n ^ m. 


P o pH 

2 ) 

<y e R 



( 

j=p 

) 


Notice that 0 ^ S(p,n), and S(p,n) is closed under multiplica- 
tion by non-zero scalars. Examples of some S(p,n) for m = 3 are 


shown : 


z 



z 



complement of this 
closed cone 
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Theorem 3 : If A is symmetric, P is a canonizing matrix for A, and 
p and n are the index and co-index, respectively, of A, then 

|x e R™ I x'^Ax > 0 = P^S(p,n)^. 

T 

Proof ; Let D be the canonical form of A, so that P AP = D = 

diag ■ ■ ,1,-1, * • ' ,-1,0,' * * ,o| with p I's and n -I's . Let x e R™ 

and let y = P ^x. Then the following statements are equivalent: 

T 

X Ax > 0 . 

T -IT -1 
X (P ) DP X > 0. 

(P“^x)'^D(P~^x) > 0. 

T 

y Dy > 0. 

2 _ 

j 

j 



p+n 

t, 

"+1 


> 0 . 


y e S(p,n) . 

X e P(S(p,n)) . 

QED. 

Corollary 5 : If W = |a^,A 2 , • • • ,A^ | is a set of symmetric matrices, 

P . is a canonizing matrix for A., p. and n. are the index and co-index 
of A^, respectively, for 1 ^ i s N, then W iS a positive definite set 
if and only if 

Pi(s(Pi,n.)) = R*" - |o} . 


Proof : This follows immediately from Definition 6 and the theorem. 

It will be recalled that a matrix A is positive definite if and 
only if every principal submatrix of A is positive definite. This 
fact generalizes to positive definite sets as follows. 
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Definition 11 : If 1 ^ t ^ m, denotes the collection of all Increas- 
ing t-termed sequences selected f rom | 1,2, • • • ,m|. If a 
In and C Is an m x m matrix, C[a] denotes the principal t x t sub- 
matrix of C lying In the aj^,a 2 »- • • ,a^ rows and a ^,02 > * • * columns of 
C. That Is, 

Cla] . . = C for l<l<t, l<j<t. 

■■ij 

Theorem 4 : Let W = |a^ ,A 2 , * ‘ * ,A^| be a set of symmetric m x m matrices. 
Then W Is a positive definite set If and only If, for each t, 1 < t < m, 
and each a c G^, = •|A^[a] , A 2 [a], A^[a]| Is a positive definite 

set . 

Proof: Suppose W is a positive definite set. Let 1 < t ^ m and a e G^., 
and let x e R^, x 0. Define y e by 


|x^ if k = ttj , j = 1,2,* •• ,t 

(O if k |a^,a2,* • • ,a^|. 


There is an 1 e |i,2,'*’,n| such that y"^A^y > 0. But then 

T t t 

X A. [a]x =5" 51 A. [a] x x = 

1 ^1 q=i 1 pq P q 


i t m m 

=1 c. >'a ■ f=l ■ 


p q p q 


The 


y A.y > 0. Thus W is a positive definite set. 

■^1 a 

reverse implication is seen by choosing a e G^, = 1 for 1 ^ 1 ^ m, 


so that W = W. 
a 


QED. 


Finally, we relate positive definiteness of a set to the number 
of positive eigenvalues of the matrices in the set. 
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I I II 


Lemma 4 ; If subspaces of R™, then 

dim 

1^1 ~ 


m. 


Proof: Let d. = dim V. and d_ = dim 
i 1 0 


(fii’i) • 


For 1 < i < N, there 


Is an (m - d.) x m matrix A. such that V. is the null space of A.. 

1 11 1 

Let A be the block matrix given by 


A = 


A is 




that is A is 


is INm - ^ For x e R^, 

Ax = 0 4^A^x = 0 for 1 < i ^ N ^x e . Thus d^ is the nullity 

of A. Now Nm - 5Zd. fe rank of A = m-d„, so 
i=l 1 0’ 

N 

d„ ^ + m - Nm . 

0 1=1 1 


That is , dim 


QED. 


Theorem 5 : If W = |Aj^,A 2 > ‘ * ,A^ | is a positive definite set of m x 
symmetric matrices and p^ = index of A^ for 1 < i < N, then 


m 


N 

JllPi > m . 

1=1 


Proof: For 1 < i < N, let P^ be a canonizing matrix for A^ and let n^ 
be the co-index of A^ . Since W is a positive definite set, Corollary 5 


gives 


Thus 


5^Pi(sCp^n.)J = r"* - {o}. 

<pi'”i>)° - hi 


c in 

, where Q denotes the complement of Q in R . 
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Now 
and P 


i 
( 6 ) 


, for each 1, S(p^,n^) c S(p^,0), so 
(s(p^,0)j ^ c P^^S(p^,n.)j^ . 

. .Thus ^ P.(s(p^,0))'= c n Pi(sCp.,np)‘= = 0 . 

i=l ' ' 1=1 ' ' 

But S(p^,0)^ ~ {^ ^ j ~ ^2 ~ ' *' ~ ^ is a subspace of 

R*" of dimension m-p^ . Since P^ is non-singular, 

Pi|s(pi,0)j ^ = P^^S(p^,0)^| is a subspace of R™ of dimens 

Thus by C6) and Lemma 4 we have 

0 = dim ^_^^P.(s(p.,0)^''j> 

( ^ (m-p . )) - (N-1) m=m-^p. ,or ^p. >m 
\i^ i=l ^ i=l ^ 


ion m-p . . 
1 


QED. 

Corollary 6 : If W = | A^, A^, • • • ,A^| is a positive definite set of 
symmetric m x m matrices, then the total number of positive roots 
of A^, A 2 , •••, Aj^, counting multiplicities, is not less than m. 
Proof: This follows. from the theorem and the fact that the index of 


A^ is the number of positive eigenvalues of A^ , counting multiplicities. 

QED. 
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VI. Conclusions. 


In an attempt to secure sufficient conditions on the feedback 

matrices F and the decision function d that will ensure stability, 

s . 

X 

we have come only this far: 

If the F are selected so that the set 
s^ 

K = <C. + D.F 
lx X s 


1 < i < 


n} 


(see section 3) 


has a positive definite Q-symmetric set for some positive definite 
matrix Q, then a decision function exists which will produce exponen- 
tial convergence (see Corollary 1) . One such decision function would 
be constructed along the lines of the algorithm given in example 1, 
section 4. 

We have also shown that if the feedback matrices are selected 
so that the resulting set is convergent, then t;he resulting set pro- 
duces asymptotic stability, and so Lyapunov stability. For conver- 


gence implies that the set 


1 < i < Nf is pre-contractlve 


<C. + D.F^ 

\ X X Si 

(Theorem 1 ') , which implies asymptotic stability (Corollary 3). 

Some of the questions which Immediately suggest themselves are: 
(in the notation of the preceding sections) 

1. Under what conditions on A, B, and- S can the F be selected to 
achieve for some positive definite matrix Q, a positive definite 
Q-symmetric set, and how should the F^ be determined? 

2. Can we find more and better characterizations of positive definite 
sets of matrices than those given in section 5? It would seem 
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that there -might be some conditions on the canonizing matrices 
of section 5 which would produce positive definiteness , 
using Corollary 5. 

3. Can we find some useful conditions relating to pre-contractive- 
ness? If K Is pre-contractive, then for some positive integer 
M, the set of all products of M or fewer members of K is con- 
tractive, so perhaps the material in section 5 will be useful 
here, too. 

4. Related to 3 Is the problem of taking a pre-contractive set K 
and forming a contractive set with finite products of members 
of K, Such a process might lead to the development of decision 
functions producing stability. 

5. In the formation of a contractive set from a pre-contractive one, 
as mentioned in 4, we may want to minimize the number of matrices 
in the resulting contractive set. Alternatively, we may want to 
minimize the number of matrices from the pre-contractive set 
which are used in forming the contractive one. It would be use- 
ful to have conditions on contractive and pre-contractive sets 
which guarantee that they are "minimal"; that is, that if any 
matrix is removed from the set, the set no longer has the desired 
property . 

. + D.F 
1 X s . 

X 

contractive, the construction of a decision function depending 
only on Xj^ is possible. In general, the decision function should 
depend only on x^ and If K is pre-contractive, we know that 

any initial x^^ can be steered to zero, but we do not yet know 


6. We have observed that if the set R = <C 
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whether there necessarily exists a decision function depending only 
on Xj^ and that can do It. If not, are there conditions on a 

pre-contractlve K which guarantee the existence of the decision func- 
tion without Implying contractiveness? 
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